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*  MOTION  CF  TRIM  BODIES  AT  HIGH  SUPERSONIC  VELOCITIES1^. 

[_  Following  is  a  translation  of  an  article  by 
G.  K* •  Bntn-Zelikovich,  A.  I..  Bunimovich  and 
M.  F<  Kikhoylove  in  |he  News  ..  of  the  Acad¬ 
emy  of  Sciences  XywwA  |  Division  of  Technical 
Sciences ,  Moscow  No*  1,  i960,  pages  33*-4cC7 
1  *Izvestiya  akademii  Nauk  SSR 

Recently,  along  with  the  theoretical  investigations  of  the 
ambient  gas  flow  past  bodies  at  high  velocities  ,  high-power  exper¬ 
imentation  facilities  have  been  built  which  make  it  possible  to 
study  the  problem  of  gas  flow  at  speeds  approve!  ng  the  velocity  of 
sound  and  considerably  exceeding  it. 

The  progress  of  the  experiment  has  attracted  the  attention  of 
the  researchers  to  the  problem ’ of  the  establishment  of  similarity  cri¬ 
teria  in  the  ambient  gas  flow  at  high  velocities  past  bodies  which  do 
net  possess  the  property  of  geometries 1  similarity* 

Ambient  go s  flow  past  thin  bodies  is  investigated  in  Kerman’s 
work  LVJ  (velocities  approaching  the  velocity  of  sound)  and  Teien’s 
L  2J  (high  supersonic  velocities),  the  flow  being  assumed  to  be 
plane  or  axially  symmetrical,  potential,  end  isen tropic,  The  consid¬ 
erations  of  these  authors  are  not  rigorous  because  of  the  assumptions 
concerning  potentiality  and  ieentropy  of 
the  flow  which  clearly  do  not  correspond 
to  the  physical  properties  cf  the  flow. 

lu  the  present-  article,  Teien’s 
results  are  generalized  for  the  case  of 
a  three-dimensional  motion  in  the  pres¬ 
ence  of  shock  waves  and  vortexes,  arid  it 
is  shown  that  the  problem  on  the  steady 

Figure  1, 

1)  The  work  was'  published  for  small  circulation  in  the  symposium 
"Theoretical  Hydromechanics" ,  No.  4,  1949*  1°  a  footnote  to  the  q  _ 

work,  it  is  pointed  out  that  an  authorized  report  of  the  papers  £ 
was  presented  in  March  1948  at.  the  seminar  on  hydromechanics  at 
Moscow  State  University*  The  work  is  being  reprinted  ithout  changes. 
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ambient  gas  flow  past  a  thin  body  at  a  high  supersonic  velocity  can 
be  reduced  approximately  to  the  problem  on  unsteady  gas  motion  in 
specs  with  the  number  of  dimensions  being  less  by  unity.  Comparison 
of  tbs  results  obtained  and  the  available  exact  solutions  determine 
the  limits  of  the  applicability  of  the  approximation  method. 

1.  Let  us  consider  the  motion  of  &  thin  body  in  the  direction 
of  the  negative  axis  x  at  &  constant  velocity  V  which  considerably 
exceeds  the  speed  of  sound  (Figure  1)* 

Equations  of  the  absolute  motion  of  gas  in  projections  on  the 
♦fixed  axes  of  coordinates  x*  ,  y* ,  z'  have  the  following  form. 

•Euler's  equations 
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Continuity  equation 


0p  ,  d  (pi/)  d  (?v)  i  dtyw)  _  q 
dt  '  dx'  1  du'  at' 


(1.2) 


Adiabatic  equation 

£  JL  +  t,±JL+w»JL 

dt  P*  a*  p*  dy>  p*  dz , 
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(1.3) 


Here,  p  is  pressure,  P  ~  density,  u,  v,  w  -  velocity  compo¬ 
nents  of  gas  along  axes  xi ,  y’* ,  z* . 

Let  us  go  over  to  the  system  of  coordinates  x,  y,  z  connected 
rigidly  with  the  moving  body,  after  having  made  use  of  the  transform¬ 
ation.  -  * 


**  =  *  —  w.  y *  s*  *»  i  (M) 


Since  the  motion  of  gas  relative  to  the  body  is  steady,  then 
for  any  function  F(x,  y,  z)  the  following  relation  is  valid 

»F  (*,  y.  *)  y  __  »F  (*' 4 Vt<  »*.  »’) 

Let  b  be  the  chord  of  the  body  and  $  -  some  linear  dimension 
characterizing  the  cross  section  (diameter  of  the  arose  section  in  the 
case  of  an  axially  symmetrical  body;  the  median  in  the  case  of  a  plane 
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'contour)*  Introducing  dimensionless  coordinates  vjt  £  according  ' 
to  formulas 


*  =  y  ~  iv4,  z—K  (16) 


end  making  use  of  the  relation  (1*5) »  we  will  obtain  from  the  equations 
(i.i)  -  (1.3) 
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If  the  form  of  the  body  is  defined  by  the  equation 

g(*/t>>  y/&,  2/5)  “g(£,TQ, C)  =0  (3.8) 

then  the  embient  flow  condition  on  the  surface  of  the  body  may  be  writ¬ 
ten  in  the  form 

<*'+">!+ ’’!+4'=°  «-9> 

Fluid  is  at  rest  at  the  infinity  ahead  of  the  body 

u  —  0,  v^O,  u>  =  0  (1.30) 

pressure  ©equals  p  and  density  equals  p0. ! 

¥hen  considering  the  motion  of  thin,  bodies  et  high  velocities, 
ve  shall  be  neglecting  terms  of  the  form  (&/ b)  u  (du  { dt),  f  b)  dp  /  d£ 

(which  contain  the  small  multiplier  (»/*)'  when  comparison, 
with  the  rest  of  the  terms  which  are  part  of  the  equations  (1.1) 

(1.2),  (1.3).  it  the  same  time ,  terms  containing  multiplier  V S/&,  are 
retained  since  at  high  velocities  V  of  the  body,  the  magnitude  of 
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V  f  b\  is  finite*. 

'  Belating  all  velocities  to  the  speed  of  sound  in  a  fluid  at 


rest,  a0,  for  which 


a°  “*P« 


Cl. 115 


(*  is  the  ratio  of  thermal  capacities),  denoting  the  dimensionless 
magnitudes  of  velocity,  density,  and  pressure  by 
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end  introducing  the  dimensionless  parameter 


we  will  obtain  from  equations  (1*?) 


(*.-  5) 


(1.12) 


(i  J3) 


Figure  2 


The.  ambient  flow  condition  (1.9)  will  be  rewritten  in  the  fol¬ 
lowing  form  after  neglecting  the  magnitude  of-u*  upon  comparison  with  V 


(US) 


Conditions  at  the  infinity  ahead  of  the  body  will  take  on  the 

form 


B*  =  0,  V*  as  0,  St?"  —  0  (1.16) 

2.  During  the  motion  of  a  body  at  high  supersonic  velocity  in 
a  ges  at  rest,  a  shock  wave  originates  which  deflects  little  from  the 
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surface  of  the  body.  If  the  normal  component  of  the  velocity  of  the  ' 
shock  wave  propagation  is  denoted  by  c  and  the  angle  between  the  P 
direction  of  the  tangent  (.Figure  2)  to  the  shock  wave  and  axis  x  is 
designated  by  £,  then  in  view  of  the  smallness  of  angle  p  we  will 
obtain 

c  =  V  sinpzzVp  (2.1) 

If  v  and  vn  are  the  components  of  the  velocity  of  gas 
along  the  tangent  end  the  normal  to  the  direction  of  the  shock  wave* 
then  rf  —  u cos p  —  w* sis pt  orbecan.se  £  is  small, 


vt  a? » 

and  V  may  be  regarded  as  the  component  of  the  velocity  of  gas  in  the 
direction  perpendicular  to  axis  x. 

Conditions  on  the  shock  wave  may  be  represented  in  the  form 


vt  -  0, 


Pi  =  P« 
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X  4-  1  *»* 


x  —  I  *$*' 
2x  c*. 


„  -  (x-M)y(x- 1) 

Pi~  Po  “i  +  2 {a<?J c*)/(x-1) 


(2,2) 


where  pp  and  p  ^  are  density  and  pressure  beyond  the  shock  wave. 

Going  to  coordinates  g,  «,  we  will  obtain  from  condition 

(2.1) 


c~V~*«  (2.3) 

where  is  magnitude  determined  in  the  process  of  the  solution  of  the 
problem  (the  tangent  of  the  angle  between  the  direction  of  the  tangent 
to  the  shock  wave  end  axis  I  in  the  plane  6,  i»j). 

Introducing  parameter  K  into  equations  (2.2)  and  denoting 
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vs  will  obtain  conditions  on  the  shock  wave  in  dimensionless  form 

,,•  =  0.  (2.5) 
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It  is  clear  from  the  equations  and.  boundary  conditions 
obtained  that  the  only  parameters  of  the  problem  are  the  magnitudes 
of  K  and  a  .  Therefore,  it  is  obvious  that  for  pressure,  deusi ty r 
and  the  components  of  the  velocity  of  gas  at  ant  point  of  the  flo*’. 


N  (2.6) 

«o  . 


3«  Tt  ia  «aay  to  show  that  in  such  an  approximate  statement 
the  problem  on  the  steady  motion  of  a  thin  body  at  a  high  supersonic 
velocity  coincides  with  the  problem  on  unsteady  motion  in  space  the 
number  of  dimensions  of  which  is  smaller  by  one , 1 ) 

By  the  substitution  of  the  variables  in  the  equations  (1,1k) 


cu:e  may  write 
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y  =  6lJ,  ■■■  Z  m,  SC 


(3.1) 


we  will  obtain  equations  of  the  non-steady -state ,  two-dimensional 
problem  in  fixed  plane  nt  perpendicular  to  the  velocity  of  the  body 
motion  (Figure  3)* 


dv  .  dv  ,  dv 

- S.  - j-  u;  —  : 

&t  ‘  dy  '  dz 

dv>  ,  dw  ,  dw 

+  vTy  +  W» 7 


dt 


dp 

at 


d  (pf) 


L 

dt' 


-f-  v 


dy 
±JL 

dy  p* 

du 

dt 


+ 


dz 
d  (pw) 
dz 


L^l 

P  dy 
dp 
P  dz 

0 


-Y  to 

=  o 


d  P 
dz  D* 


(3.2) 


0 


Condition  (1,15)  goes  into  the  equation  of  the  motion  of  the 
boundary  which  may  be  regarded  as  an  unusual  piston.  Indeed,  since 


d&  __dg  b_  d£  _  dg  . 

flf  dtV  ’  drj  dy  ’  3%  8z 


(3.2a) 


1)  The  feasibility  of  reducing  the  problem  on  the  three-dimensional 
steady  motion  (in  this  approximate  statement)  to_ the^two-dimensicnal 
problem  on  unsteady  motion  is  pointed  by  Keyes  L>  J- 
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v  **  0,  u;  =  0  ii | mi  t  s=  0  (3.4) 

'  The  'change  in  £  under  the  condition  that  K  remain  constant 
is  simply  a  change  in  the  scale  in  the  non-steady-state  problem. 

Thus,  the  problem  on  the  motion  of  a  thin  body  at  high  super¬ 
sonic  velocity  is  three-dimensional  space  corresponds  in  the  approxi¬ 
mate  statement  to  the  plane  non-steady-state  problem  of  the  expansion 
of  * 'cylindrical  piston. 

The  plane  problem  on  the  motion  of  a  thin  body  at  a  high  super¬ 
sonic.  velocity  corresponds  to  one-dimensional  problem  on  the  unsteady 
motion  of  a  piston  in  the  presence  of  the  shock  wavs. 

it.  In  order  to  determine  the  limits  of  the  applicability  of  . 
the  similarity  criteria  obtained,  let  us  compare  the  approximate  solu¬ 
tion  with  the  available  exact  solutions. 

As  the  first  example ,  let  us  consider 
the  problem  on  the' motion  of  a  wedge  having 
angle  2€>  in  the  direction  of  axis  -x  at 
supersonic  velocity  (Figure  k) « 

As  is  known,  a  straight -line  shock 
wave  forms  in  the  flow,  this  shock  wave 
originating  ,at  the  point  of  the  wedge  and 
forming  angle  f,  with  axis  x. 

Figure  A 

It  can  b©  derived  from  the  conditions  on  the  shock  wave  (2,2) 

that 


„  o  ?L<*  -  tt  AV  sir  &.+  2}  tg  ft 

w  —  P  —  »«cig  — <*qri}  Mft  sif}5  p 

Pi  =  [2xAf0*  sin*.p  —  (*  —  i)J 


whence  it  is  possible  to  find  the  relation  P^/ p9  as  the  function  of 
the  number  MQ  and  the  central  angle  Q  of  the  cone.  In  the  ease 
being  considered*  parameter  K  equals 


tge  (4.2) 
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Figure  5 


Graphs  of  the  relation  of  p-J pQ 

to  K(x.ziJ>)  h era  been  constructed 
(Figure  5)  on  the  basis  of  the  com¬ 
putation  performed  for  numbers  MQZ 
Z  1*5*  2.-5 *  10  and  variations  of 
angle  H  from  0  to  20O, 

The  approximate  solution  of  the 
problem  on  wedge  motion  at  high 
supersonic  velocities  amounts  to  the* 
solution  of  OD.e-dio*ensiorial  problem 
on  piston  motion  at  a  constant  speed. 
The  speed  of  piston  motion  is 
obviously  equal  to  the  vertical  ve- 
loci ty  of  the  wedge ,  he. 

v  =  F  tg  0  (4.3) 

and  consequently, 

ge-2-f  (4.4) 


as  is  known,,  in  piston  motion  the  pressure  and  velocity  on 
the  piston  ©qua.,  the  pressure  and  velocity  beyond  the  shock  wave. 
Eliminating  the  magnitude  at)/e  from  the  second  and  third 
equations  of  (2,  .2)  and  making  use  of  (4*4)  *  we  will  obtain 
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The  curve  corresponding  to  (4,5)  is  represented  in  Figure  5 
by  the  broken  line* 

It  is  evident  from  Figure  5  that  for  large  numbers  for  M0 
and  smell  central  angles  of  the  wedge,  the  approximate  solution 
coincides  very  well  with  the  exact  solution* 

For  comparison  there  are  plotted  in  Figure  6  broken-line 
curves  which  correspond  to  the  approximate  computation  of  Pq/p&  by 
Bonov's  formula  £~k  J  if  one  confines  oneself  in  it  to  the  terms 
of  the  second  order  relative  to  0.  Doaov«s  formula  has  the  form 
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where  K  =  2Mo0is  taken  approximately 


Figure  6.  Figure  7* 

&b  the  second  example,  let  us  consider  the  motion  of  a  right  circu¬ 
lar  cone  in  the  direction  of  its  axis  at  a  supersonic  velocity. 

The  exact  solution  of  this  problem  was  given  by  Busemann 
Z*3  J  and  the  detailed  computations  were  performed  by  Taylor  and 

Mac  Coll  [_  6 

¥e  showed  above  that  this  problem  corresponds  approximately 
to  the  problem  on  the  expansion  of  a  circular  cylindrical  piston 
at  a  steady  velocity.  And  the  solution  of  the  problem  on  the  expan¬ 
sion  of  a  cylindrical  piston  vas  given  L*  I*  Sedov  /  ?  J. 

In  Figure  7  are  given  the  curves  of  the  relation  between  the 
ratio  of  the  pressure  on  the  cone  (cylindrical  piston)  to  pres¬ 

sure  p0  in  gas  at  rest and  the  number  for  K  sz  (V  ] a0)  2  tg  0  —  2  vjao 
(because as  in  the  case  of  the  wedge,  S/b*2tg0,  and  the  velocity 
of  piston  motion  is  the  vertical  cone  velocity  equalling  Vtg  0). 

The  broken  line  represents  the  result  of  numerical  integration  of 
equations  derived  in  Sedov's  work.  Continuous  lines  show  the  depend 

enee  of  Pk/p0  on  K  with  different  numbers  for  V/aQ,  this  rela¬ 
tion  having  been  computed  on  the  basis  of  the  exact  solution  by 
Taylor  and  Mac  Coll*  Dot-and-dash  lines  show  the  dependence  of 
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2%/$0  on  K  with  different  numbers  for  ^o/»0  *  this  relation 

having,  been  obtained  by  Basemana*  a  formula  £  8  J- 

It  ia  evident  from  Figure  f  that  the  difference  in  the  values 
produced  by  the  exact  end  approximate  solutions  decreases  with,  the 
increase  in  the  number  for  M  .  For  Mei3*  the  result  obtained  from 

the  approximate  solution  may  be  considered  satisfactory  when  K  < 1 
(an  error  of  less  than  $%) *  i.e*  for  the  central  angles  of  the  cone 
,  2  ft  <  18°.  For  Mcr5*  the  result  is  satisfactory  up  to  KsS*  i*e. 
for  cones  with  central  angles  2  ,4)  <  22°.  Generally*,  however ,  pres¬ 
sures  computed  by  approximation  method  are.  greater  in  magnitude  than 
pressures  computed  by  the  exact  method.  Conversely*  pressures  comput¬ 
ed  by  linearized  method  are  smaller  in  magnitude  than  pressures  com¬ 
puted  by  the  exact  method*  the  error  increasing  with  the  increase  in 
the  number  for  KQ, 

5»  A«  has  been  shown*  with  the  assumptions  made  above  concern¬ 
ing  the  form  of  the  bodies  and  the  flow  velocity*  all  dimensionless 
dynamic  elements  of  motion  depend  only  on  the  dimensionless  parameter 
K»  Making  us©  of  this,  let  ua  find  now  the  functional  relationship 
for  the  coefficients  of  head-on  drag  cx  and  lift  fores  cv„ 

ife'i ng  of  Infinite  Span.  Let  us  define  the  equation  of  contour 
in  the  form 

y  =  $//(!)  =  $#(*.)  (5.1) 


where  R(  )  is  the  distribution  function  of  the  wing  thicknesses. 
Than,;  making  use  of  (2.6),  we  will  obtain  for  the  aggregate  drag  and 
lift  force 

X  «  ^  p  g  dx  *  8  ^pH '  ($)  $  =  *p$i  (&)  (5.2) 


K  —  ^  pdx  — - 

'  -  1 


(5.3) 


Expressions  for  the  coefficients  of  drag  and  lift  force  have 
the  form 
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(5.5) 
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For  bodies  similar  ia  the  sense  indicated*  Akkeret’s  linear¬ 


ized  theory  gives  the  coefficients  of  cx  and  cy  ia  the  form 


(6/6)* 


(6  fb) 

■»  W  -  A~-r-:r — ~~~T~ 


(5.6)  [ 


For  large  values  of  the  number  for  Mc*  these  expressions 
assume  the  form 


(6/6)*  __  Jf* 


AV  * 


,  _<*/*>  _  *_ 
Af,  A4* 


which  agrees  with  the  equations  (5*4)  end  (5*5)  although  these  .equa¬ 
tions  are  considerably  more  general. 

Axially  Symmetrical  Body.  Let  the  equation  of  the  body  be  de¬ 
fined  in  the  form 

r  *  Ui  (5)  «  8 h  (jj  (5.8) 

where  is  the  maximum  diameter  of  the  cross  section;  m-  the 
distribution  function  of  thicknesses,.  Then,  we  will  obtain  for  the 
total  drag 

X  -  U2.r*  =  2nS V„  \  £*<*)*' ©.ft  =  2.SWWK)  (5.9) 


The  coefficient  of  drag  related  to  the  area  of  the  maximum 
cross  section  is 


,  2 vWp&km  r*(K) 

c* "  (Vdi^ra^ 


(5.10) 


The  similarity  laws  obtained,  show  that  for  bodies  with  the 
same  distribution  of  thicknesses  at  the  angles  of  attack  proportional 
to  the  relation  8  /  b,  amounts  cxM03  and  CyM02  (wing  of  infinite 

span),  (axially  symmetrical  body)  will  be  functions  csf  one 

parameter  K $i0bjb. 

Consequently,  having  experimental  results  of  blowing  or  of 
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flight  tests  for  some  profile  with  different  numbers  for  Mc*  it  is 

possible  to  recompute  the  obtained  results  for  a  series  of  prefiles 
with  the  same  distribution  of  thicknesses.  Conversely*  after  blowing 
upon  a  series  of  profiles*,  similar  in  the  sense  indicated,  with  the 
seme  number  for  M0*  it  is  possible  to  recompute  the  results  for  each 
of  the  profiles  of  the  series  with  different  numbers  for  Mc« 

Formulas  derived  in  the  present  article  for  the  magnitudes  of 
cx  and  Cy  after  taking  into  account  the  presence  of  shock  waves 

and  vortexes  in  the  flow,  coincide  with  Teien’s  results  obtained 
with  the  assumption  of  potentiality  and  isentropy  of  the  flow. 

Submitted  26  August  1959 
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